Implementation of Clifford gates in the Ising-anyon topological quantum computer 
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We give a general proof for the existence and realizability of Clifford gates in the Ising topological quantum 
computer. We show that all quantum gates that can be implemented by braiding of Ising anyons are Clifford 
gates. We find that the braiding gates for two qubits exhaust the entire two-qubit Clifford group. Analyzing the 
structure of the Clifford group for « > 3 qubits we prove that the the image of the braid group is a non-trivial 
subgroup of the Clifford group so that not all Clifford gates could be implemented by braiding in the Ising 
topological quantum computation scheme. We also point out which Clifford gates cannot in general be realized 
by braiding. 
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I. INTRODUCTION 

Quantum Computation (QC)i development encountered 
tremendous difficulties in storing and manipulating quantum 
information in real physical systems because of the over- 
whelming decoherence and noise. Topological Quantum 
Computation (TQC)^ is a branch of QC in which both infor- 
mation storage and processing are protected by the topolog- 
ical nature of the quantum computer In the TQC approach 
the quantum information is encoded in non-local topologi- 
cal degrees of freedom and is therefore inaccessible to noise 
and decoherence which are mainly due to local interactions. 
Moreover, the quantum gates are implemented by non-trivial 
topological operations which are once again protected against 
decoherence. The purpose of this topological protection of 
qubits and quantum gates is to improve the quantum compu- 
tation hardware to such an extent that quantum information 
processing become more feasible in real physical systems. 
The feasibility of this approach depends on the detailed noise 
structure^. In this paper we work entirely in the noiseless case, 
and ask about the possibility of realizing certain key gates in 
a specific scheme. 

One of the most promising TQC schemes^ employs the an- 
ticipated non-Abelian braid statistics of the lowest energy ex- 
citations, called Ising anyons, of the fractional quantum Hall 
state at filUng factor v = 5/2, which is believed to belong to 
the universality class of the Moore-Read (MR) Pfaffian stata^. 
Due to the topological degeneracy in the two-dimensional crit- 
ical Ising model, representing the neutral degrees of freedom 
in the MR state, it becomes possible to realize «-qubits by 
2n + 2 Ising anyons: the states of 2n + 2 Ising anyons are 
represented by conformal field theory (CFT) coiTelation func- 
tions (more precisely, chiral CFT blocks) which happen to be- 
long to one of the two inequivalent spinor irreducible repre- 
sentations (iRs)6j7j8j9ji0jii_ q£ j-jjg covering group Spin(2« + 2) 
of the rotation group SO(2n + 2). There are two inequiva- 
lent spinor lR aV.io,ii so(2n + 2), of dimension 2", which 
differ by their total fermion parit y Despite being math- 
ematically inequivalent they appear to be equivalent from the 
computational point of viewii, i.e., the set of matrices that 
could be obtained by braiding Ising anyons in both represen- 



tations is the same. 

Clifford-group gates, that are defined as those unitary op- 
erations that preserve the Pauli group, play a central role in 
quantum information theory. Although the Clifford group is 
not sufficient for universal QC and its computational power 
cannot exceed classical computers'''?-'- , experimental real- 
ization of the Clifford group with sufficient scalability would 
be a cornerstone of QC. Indeed, the Clifford group is of sig- 
nificance for quantum error correctiorU^ and allows the gener- 
ation of entangled states, like GHZ or Cluster states, the latter 
being a prerequisite for universal quantum computation in the 
measurement-based schemai^. 

Using the explicit representation of the braid-group gen- 
erators for the exchanges of Ising anyons one of us has 
constructed^ the entire Clifford groups for 1 and 2 qubits in 
terms of braid generators for S§n and ^g, respectively. How- 
ever, this approach encountered serious difficulties for em- 
bedding some Clifford gates in systems with three or more 
qubits^: the topological entanglement between distant Ising 
anyons induces additional phases when exchanging anyon 
pairs which are in the state i.e., when exchanging pairs 
sharing Majorana fermions. Therefore, it was possible to con- 
struct only a part of the Clifford group for 3 qubits. In this pa- 
per we shall address the question whether all Clifford-group 
gates could be realized by braiding of Ising anyons in the Ising 
TQC scheme or not. Because, as we shall prove below, the 
«-qubit Pauh group coincides with the monodromy subgroup 
representation for 2« + 2 Ising anyons and because in gen- 
eral the monodromy group is a normal subgroup of the braid 
group, it naturally follows that the Ising-model braid group 
representation is a subgroup of the Clifford group. In other 
words, all quantum gates that could be implemented by braid- 
ing of Ising anyons are actually Clifford gates. Unfortunately, 
it also appears that not all Clifford gates could be realized by 
braiding for three or more Ising qubits and we shall try to ex- 
plain why. The Clifford gates that cannot be implemented by 
braiding are typically the embeddings of the two-qubit SWAP 
gate into larger systems. 

The rest of this paper is organized as follows: in Sect. 
we describe how anyonic states of matter could be labeled 
by fusion paths in Bratteli diagrams and how this could be 
used to determine the dimension of the computational space. 
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In Sect. Hn] we summarize the explicit representation of the 
elementary Ising-model exchange matrices as proposed by 
Nayak and Wilczek— . In Sect. |lV]we explain the symplectic 
description of the Clifford group for « Ising qubits and esti- 
mate the order of the Clifford group in order to compare it 
with the order of the braid-group representation. In Sect. [V] 
we analyze the relation between the Pauli group and the mon- 
odromy subgroup of the braid group and explain why all quan- 
tum gates that could be realized by braiding of Ising anyons 
are in fact Clifford gates. In Sect. IVIII we give the explicit 
braid construction of the two-qubit SWAP gate which allows 
to construct the entire two-qubit Clifford group by braiding. 
Some important technical details are collected in several ap- 
pendices. 



n. FUSION PATHS: LABELING THE ANYONIC STATES 
OF MATTER 

The anyonic states of matter (should they really exist in Na- 
ture) differ from the ordinary fermionic and bosonic states in 
that we need to give additional non-local information in or- 
der to specify the quantum state. In TQC this extra infor- 
mation is expressed in terms of topological quantum num- 
bers that are eventually used to encode quantum information 
non-locally in order to gain topological protection. An impor- 
tant ingredient of any TQC scheme is the necessary degen- 
eracy of ground states of the multi-anyon system in presence 
of trapping potentialsi^sii (the potentials that keep our com- 
putational anyons at fixed positions). This is equivalent to 
a degeneracy of the multi-anyon states, considered as exci- 
tations corresponding to having a number of anyons at fixed 
positions in the plane, over the ground state (this time with- 
out trapping potentials In the CFT language, that we will 
use to characterize the anyonic states as CFT correlators, the 
second point of view is more appropriate and we shall speak 
about the degeneracy of the CFT correlation functions corre- 
sponding to 2n + 2 Ising anyons at fixed positions. In more 
detail, the CFT correlation functions of 2n + 2 Ising anyons at 
fixed positions inside a Pfaffian droplet^ span a Hilbert space 
of dimension 2" and are therefore appropriate for representing 
n Ising qubits. The non-local internal quantum numbers, for 
non-Abelian anyons, that would allow us to distinguish be- 
tween the different states in this computational space are the 
so called fusion channels. In order to make this notion more 
transparent let us try to explain the fusion rules of the Ising 
model: the most relevant quasiparticle excitations in the low- 
temperature, low-energy regime, are described by the chiral 
Ising spin field a with dimension 1/16. If two such quasipar- 
ticles, having non-trivial topological properties, are fused, i.e., 
taken to the same point in the coordinate plane the composite 
object would look (from far away) like another topological ex- 
citation and this is symbolically expressed by the fusion rule 



ax a 



(1) 



prepared experiments some pairs of a would behave collec- 
tively as the vacuum (this is the I in Eq. ([T])) and this is called 
the vacuum channel, while some others would look like the 
Majorana fermion (i.e., the xjr in Eq. ([T])) and this is called the 
Majorana channel. In other words, the combined quasiparti- 
cle can be considered as a mixed state of the vacuum and a 
Majorana fermion. Localized particle-like collective excita- 
tions, such as the field a in Eq. ([T]i, which have more than 
one available fusion channels are called non-Abelian anyons. 
There are in general superselection rules which forbid creation 
of coherent superpositions of anyons, such as the vacuum and 
the Majorana fermion in the Ising model, belonging to differ- 
ent superselection sectors and this has to be taken into account 
when using non-Abelian anyons for TQC. 

The important point is that the topological properties of the 
anyon pairs are persistent (under some reasonable assump- 
tions), i.e, if two quasiparticles, which are in a state character- 
ized by a definite fusion channel, are pulled away the pair still 
possesses the properties of their corresponding fusion channel 
and, e.g., if they are fused again after some time, they will pro- 
duce the same result as that determined by the original fusion 
channel. 

When we have many Ising anyons, which we will assume 
to be ordered on a line, we could represent the string of 2n + 2 
Ising anyons (7 into n + l pairs (<7, <7)e and characterize each 
pair by its fusion channel c. Then the sequential composition 
of the fusion channels for all pairs can be described by a fu- 
sion path in the corresponding fusion diagram. The fusion di- 
agrams that we will use are known as Bratteli diagra m'^''^-'^ 
and represent graphically the possible results of fusion of 
a single basic non-Abelian anyon to an array of other non- 
Abelian anyons usually of the same type. One way to under- 
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This rule means that there are two distinct channels in the fu- 
sion process and if we consider a large number of identically 



FIG. I: Bratteli diagram for the Ising model, in the su{T) notation. 

stand this diagram is to represent the Ising model as the affine 
coseti^ iM(2)2/M(l) and use the fact that the Ising model fu- 
sion rules are the same as those for the su{2)2 Wess-Zumino- 
Witten model Then the CFT primary fields a and i// (of 
CFT dimensions 1/16 and 1 /2 respectively), could be labeled 

by the reduced Young tablues for the admissibloi^ ■s"(2)2 rep- 
resentations as shown in Fig. [Hand define (together with the 
vacuum) the 3 different superselection sectors (anyons) of the 
Ising model. In this way the fusion rules for the Ising model 
can be inferred from the tensor product decomposition of the 
su(2) representations. Each step to the right in Fig. [T] denotes 
fusing one more fundamental anyon a to the existing string 
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of (7 anyons (whose length is determined by the x coordinate) 
and the arrows represent the possible fusion channels for this 
process, which are listed along the vertical axis. Note that the 
Bratteli diagram in Fig.[T]is finite in vertical direction which 
expresses the important property of rationality: any rational 
CFT, such as the critical two-dimensional Ising model, con- 
tains a finite number of topologically inequivalent superselec- 
tion sectors (i.e., finite number of distinct anyons) which are 
"closed" under fusion. In order to build the Bratteli diagram 
in Fig. [T]we have to supplement Eq. ([T]i with the two other 
fusion rules y/x<7 = <7, Ix<7 = <T. All anyons, such as I and 
i/A, from which originates only one arrow pointing to the right 
are Abelian, while those, such as a, with more than one (two 
in this case) arrows to the right are non- Abelian. 

In order to gain some intuition about labeling the states in 
the Ising model we show in Fig. |2] two distinct three-qubit 
states realized as fusion paths in the Bratteh diagram for 8 
Ising anyons. Notice that there are two inequivalent repre- 
sentations of the multi-anyon states with 8 anyon a^i^*^''^ : be- 
sides the 8 <J fields inside the CFT correlation function there 
are also a number of Majorana fermions present in there (not 
shown explicitly in the correlation functions in Fig. |2]but see 
E q. (|4|i below) - when this number is even the state belongs 
to the positive-parity representation (in which case the 8 (7 
fields must fuse to the identity in order for the CFT correla- 
tor to be non-zero), while when it is odd the state is in the 
negative-parity representation (and the 8 <J fields must fuse to 
the Majorana fermion). The Bratteli diagrams also appear to 



0) 
c 
c 
CO 
.c 
o 
c 
o 



I 



((o-<7) (a<7U<7<7),(<7a),)^|lOO). 
((oD-y (crcr) (crcr) (aery =|011 




number 
of cTs: 



FIG. 2: (color online) Fusion paths in Bratteli diagram: red circles 
display a path corresponding to the three-qubit state |100)_, while 
the path of blue squares corresponds to the state |011) + . One hor- 
izontal step to the right represents adding one more a field and the 
arrows point to the corresponding fusion channel. The black circle 
denotes the vacuum channel I. The subscript ± of the computational 
states denotes the positive- or negative-parity representations. 

be a very useful graphical tool for computing the dimensions 
of the computational spaces for TQC with anyons. For exam- 
ple, the space dimension of Ising-model correlation functions 
with 8 anyons at fixed positions could be read off from Fig.|2] 
as follows: each new step to the right defines a single new fu- 
sion channel if the current number of a is even and two fusion 
channels if the number of as is odd. Therefore, the number of 
distinct fusion paths is just 2t , the power of two with an 
exponent equal to the number of the odd c7-steps and finally 
divided by 2 because the last step is always fixed to fuse to the 



vacuum or to the Majorana fermion depending on the parity 
of the representation we deal with. We see from Fig. |2]that the 
number of independent correlation functions of 8 fields a is 
2^^ ' = 8 which could be easily generalized to the case of 2n 
a fields where the correlation functions (for fixed positions of 
the anyons) span^ a computational space of dimension 2"^^ . 

Thus we can say that the state of a multi-anyon system can 
be ultimately characterized by its fusion path in the Bratteli 
diagram. Having specified a multi-anyon state of matter as a 
fusion path in a Bratteli diagram any TQC scheme requires 
physical processes that could initialize the multi-anyon sys- 
tem into a given n-qubit state and this has to be further sup- 
plemented by procedures for measuring the states of the indi- 
vidual qubits. For the TQC scheme based on Ising anyons 
the initialization as well as the measurement by a Fabry- 
Perot or Mach-Zehnder interferometer have been discussed 
in Refs. HielfTTL 

Now that we have got an idea about how to encode 
qubits using non-Abelian anyons we could think about exe- 
cuting quantum gates which are unitary transformations act- 
ing within the computational space of multi-anyon states la- 
beled by fusion paths. A central issue in this context is 
that the multi-anyon states, that we would like to use to en- 
code qubits, are degenerate in energy (at least approximately, 
see Refs.[lil3 for more detailed explanation), and separated 
from the rest of the excitation spectrum by a gap. This al- 
lows us to apply a version of the adiabatic theorem which 
is appropriate for the degenerate case. In simple words the 
adiabatic theorem in this case states that if the initial state is 
in the degenerate subspace that is separated by a gap from 
the rest of the spectrum and we consider the adiabatic evo- 
lution, when some of the anyons are transported along com- 
plete loops around others, the final multi-anyon state would 
be again a member of the same degenerate subspace (e.g., of 
ground states in presence of trapping potentials). Then the va- 
lidity of the adiabatic approximation would guarantee that the 
transformation of the initial state into the final one is described 
by the action of a unitary operator which includes the Berry 
phase and the explicit monodromies of the (typically multi- 
valued) multi-anyon states. It is, however, possible to choose 
a basis of CFT blocks^ in which the Berry phases are trivial 
and the entire effect of the adiabatic evolution is contained in 
the monodromies of the CFT correlation functions. Thus we 
see that it might be possible to execute quantum gates over 
our topologically protected qubits by adiabatically transport- 
ing some anyons around others and these quantum operations 
are naturally protected against noise and decoherence. 

Using the CFT-correlators representation we can describe 
the «-qubit Ising system by a CFT correlation function includ- 
ing 2n + 2 Ising spin fields a and the quantum-information 
encoding rule could be chosen to be e.g.. 



(((7(7),, • • • ((7(7),, • • • ((7(7),„ ((7(7),o)cFT 



where C{) = c\C2 - ■ ■ c„. Notice that the n qubits are encoded 
from left to right starting with the first pair of a fields and 
there is always one extra pair, which we shall take as the right- 
most one, which is inert, i.e., it contains no independent in- 
formation because its fusion channel cq is determined by the 
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product of the fermion parities of the individual qubits states. 
In other words, the physical meaning of c, is the fermion par- 
ity of the pair representing the /-th qubit, the role of the inert 
pair is to compensate the combined fermion parity of the first 
n pairs so that the correlation function is non-zero. 

Since we shall use the fermion parity of the non-Abelian 
CJ fields to encode information in the Ising qubits it is worth 
saying a few words about it. The chiral spin fields a are 
identified with the primary fields of the Ising CFT that in- 
tertwine between the vacuum sector (or the Neveu-Schwarz 
(NS) sector) and the so called Ramond (R) sector of the Ising 
model^'^'^'. The conservation of the fermion parity, combined 
with the fact that the Majorana fermion has a zero mode in 
the R-sector, implies that Ramond sector in the Ising model is 
double degenerata^Si^. This means that there must exist two 
Abelian chiral spin fields a+ and a of CFT dimension 1/16 
with fermion parity + and — respectively. However, this dou- 
ble degeneration of the R-sector is incompatible with the mod- 
ular invariance which is at the heart of the fusion rules for the 
anyonic model. In order to guarantee the modular invariance 
of the model we need to choose only one linear combinational 
of C7+ and O which is conventionally taken as 



C7 



V2 



(2) 



As a result of this (GSO) projection--' the chiral fermion par- 
ity is spontaneously broken which is also obvious in the fu- 
sion rules ([T]). Note that conservation of chiral fermion parity 
requires that the fields a± with definite parity are Abelian, 
i.e., their fusion rules must be a+ x a+ = I, (7_ x (7_ = I 
and C7+ X (7- = i//, and the non-Abelian statistics appears only 
when this symmetry is broken (which could only be sponta- 
neous because the generator of the fermion parity symmetry 
commutes with the conformal Hamiltonian). 

Although the only field which is believed to appear in the 
physical system is (|2]i and the chiral fermion parity is broken 
by the non-Abelian fusion rule ([T]) the fields a± are still con- 
venient for labeling our computational basis in terms of Ising- 
model correlation functions. In other words, even if a± are 
unphysical they could efficiently label the independent func- 
tions spanning the degenerate space of correlation functions 
because the product e i = c = ± 1 of the two indices in a pair 
(Tg, Ge^ determines the fusion channel c ~ +1 for the vacuum 
channel and c = — 1 for the Majorana one. Thus, the general 
qubit encoding scheme for Ising qubits can be represented as 

|0) ^ (<7C7)+i <^ <7+C7+ and 

|1) <^ (ccr)-! <;=^ c7+(7_. (3) 

More precisely, using for instance the CFT-description of 
Ising qubits in the positive-parity representation, and intro- 
ducing one extra pair as illustrated in Fig. |3] we could write 
the single-qubit computational basis as 

2N 

|0>+ = {a+a+ a+a+Yl\i/izj))cFT 
j=i 

IN 

= {a+a- a+a-Y\\if{zi))cFi, (4) 




FIG. 3: Single qubit configuration in terms of 4 Ising quasiholes 
coffesponding to the positive-parity representation 5+ of the braid 
group .^^4. 



where the first pair of o fields represents the qubit while the 
second one is the non-computational inert pair of Ising anyons 
compensating the fermion parity. The total topological charge 
of the 4 a fields is +1 which corresponds to the positive- 
parity spinor representation of SO(4) (cf. next section) and 
is denoted as a subscript of the computational basis states in 
Eq. (|4|. The computational basis for two Ising qubits is explic- 
itly written in Sect. lVIll and for more qubits in Refs. flolfTll 



III. ELEMENTARY EXCHANGE MATRICES: SPINOR 
APPROACH 

As we have already mentioned the n-qubit system could 
be represented by 2« + 2 Ising anyons because the degener- 
acy of the states containing 2« + 2 anyons is^ 2~^^' = 2". 
The most natural approach to determine the braid matrices 
representing the exchanges of adjacent anyons over the (de- 
generate) space of correlation functions is to use the analytic 
continuation of these multi-valued analytic functions that have 
been used to define the qubits. The results of Ref.^ suggested 
that the multi-anyon Ising wave functions span a spinor repre- 
sentation of the rotation group SO(2n-|-2) (more precisely, 
one of the inequivalent spinor representations of its cover- 
ing group Spin(2« + 2)), however, this important conjecture 
has not been completely proven in Ref. ^ It has later been 
proven rigorously in Ref. |9| for the positive-parity representa- 
tion of the braid group ^4 and in Ref. [l^ for both positive- 
and negative-parity representations of the braid group .'^2n+2 
representing the exchanges in the «-qubit states defined as 
CFT correlation functions of 2« + 2 Ising anyons. One of 
the main results of Ref. [l^ is that the «-qubit representations 

of the braid-group generators obtained by analytic 

continuation of the CFT correlation functions, are completely 
equivalent to those, denoted as 7^*."+''*'^ obtained from the 
standard Clifford-algebra realization of the generators of the 
rotation group SO(2« + 2), i.e.. 



B 



(2h+2,±) 



^(^(2h+2,±)^ ' ^(n+l,±)(-.(2n+2,±) 



Notice that this is not completely trivial since the matrices 
which establish the equivalence between the braid genera- 
tors in the wave-function representations with definite parity 
and those in the Clifford-algebra realizations with the same 
parity are special non-singular entangling matrices (i.e., ma- 
trices which cannot be written as factorized tensor products 
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of smaller ones) which are explicitly constructed in Ref. [T^ 
For example the equivalence matrices for 3 qubits in the two 
representations are C^^'^^ — diag(— 1, 1,1,1,1,1,1,-1) while 
— /diag(l, 1, 1,-1, 1,-1, — 1,-1). Thus, as proven in 
Ref. [l^ it is completely legitimate to use the Clifford-algebra 
construction of the braid-group generators for those represent- 
ing the coordinate exchanges in the CFT correlation functions. 
The advantage of using the Clifford-algebra representation of 
the braid generators comes from the possibility to express the 
Pauli gates in terms of the squares of elementary braid genera- 
tors (cf. section |V]i and then use the group-theoretic approach 
to find the stabilizing Clifford group. 

That is why we shall review below the general Clifford- 
algebra construction of the two inequivalent spinor represen- 
tations of the braid group 3^2ii which will be used later to 
implement various quantum operations over the Ising multi- 
qubit systems. Nayak and Wilczek suggested^ that the ele- 
mentary matrices representing the coordinate exchanges of 2n 
Ising quasiparticles in the Pfaffian FQH wave functions can 

be expressed in terms of the gamma matrices y- , 1 < / < 2n, 
which satisfy the anticommutation relations of the Clifford al- 
gebra 



{^W,7r}=25,„ l</,;<2„. 



(5) 



In more detail, the elementary operations for the exchange of 
the j-th and {j + l)-th quasiparticles could be expressed (in- 
cluding the correct phase) a a^i^^i^'i-^^ 

/^r=e'fexp(-|y!")r!:'0-i(^-^M:'0- 

The 2n matrices jj"^ have dimension 2" x 2" and can be de- 
fined recursively in terms of the Pauli matrices a^, = 1 , 2, 3) 
as follows^ 



(n+l) 
hn+l 

{n+l) 
hn+l 



7/ '® (^3, 
I2" ® CTi 
I2" CSi (Jl- 



1 < < 2n 



(7) 



Starting with « = as a base, where /I'' = Ci and /t'' = CJ2 
we could write the gamma matrices explicitly as follows^^f^l 



7i = (Jl (g) (J3 (8) • • • (8) (T3 
72"' = (T2 (8) 0-3 ® ■ ■ ■ ® 0-3 

V ' " V ' 

= h<»- - -® I2 g2 g3 (g) ■ ■ ■ (g) g3 

j-1 n-J 



(«) 



1211-1 ' 



)ai 

)(J2. 



(n) 

The "gamma-five" matrix 7^ could be constructed by 



(n) / .\fj (n) (n) 

7f -(-0 7i •••72„ =(^3( 



)(T3. 



Because anticommutes with all jj" it commutes with all 
Rj and therefore this representation of the braid group is re- 

(n) 

ducible. Since has only two distinct eigenvalues, ±1, it 
splits into two irreducible representations that could be ob- 
tained with the help of the two projectors 



pW=.!^,where(pW)^ 



(9) 



In other words, the generators of the two irreducible spinor 
representations of the braid group ^2n can be obtained by 
simply projecting Eq. (|6]l, i.e.. 



7?("-±) = e'fpWexp(-^7;7/+i)^r 



7r_ 
4 



(10) 



It would be useful for our analysis to give a more detailed 
expressions for the braid generators 

Rf._, = Hl(l2-/cTi^))pW, l<j<n 



e'T 



?± 

" V2 



= ^ (I4 - /cj^^') pi"', 1 <;•<„-! 



where we have used the shorthand notation CT^'''' to denote 
acting on qubit /'. A more explicit form of the matrices can be 
found in Ref. [nl 



(8) 



IV. CLIFFORD GROUP FOR n QUBITS: SYMPLECTIC 
DESCRIPTION 

Quantum operations belonging to the Clifford group are 
characterized in the Heisenberg picture as the group of unitary 
operators which stabilize the Pauli group, i.e. they map an ar- 
bitrary factorized tensor product of Pauli operators to another 
factorized tensor product of Pauli operators. This property 
has important consequences, extending a suitable description 
of the Pauli group (cf. Eq. ( ITtI i) to the Clifford group equips 
us with an efficient way of writing down such unitaries (cf 
Eq. (fTsl l) as opposed to general unitaries where the number 
of parameters typically grows exponentially with the number 
of qubits. Although, this eliminates exponential speed-up, the 
Clifford group is of great relevance for quantum computation. 
On the one hand we can study basic properties of quantum me- 
chanical systems, like entanglement generation, explicitly and 
on the other hand there are ingredients of quantum computa- 
tion, e.g. error correction'^, which are based on the Clifford 
group. Detailed descriptions of Clifford group quantum com- 
putations are given in Ref. [llll2il3[ In the following we will 
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introduce an isomorphism ( fTSl ) between a factor of the Clif- 
ford group and the symplectic group Sp2n{2), which will be 
used to further analyze the structure and computational power 
of the Ising model TQC (cf. appendices iB] and ICli . 

The Pauli group ^„ for « qubits is generated by the Pauli 
operators a, acting on each of the qubits. The elements of ^„ 
can thus be written as tensor products of the Pauli matrices 
multiplied by overall phases, which are powers of /, 



{ 



I"' a, 



a(l) 



a{j), me {0,1,2,3}}, 

(11) 

where we have defined Ob =^2- The Clifford group %, for n 
qubits is defined as the normalizer of the Pauli group, i.e., the 
group of those n-qubit unitary transformations U which map 
the Pauli group to itself under conjugation 



f„ = {u <ESu{r) I u* 



■Ud^n }. 



(12) 



Clearly, '^n is an infinite group since if U satisfies Eq. ( fT2] ) so 
does e"^U for any G M. However, the overall phases e"^ are 
in general irrelevant for quantum computation and therefore 
we could remove the center Z of the Clifford group, which is 
generated by these phase factors. The resulting factor group 

which we will call the projective Clijford group, is already a 
finite group whose order will be derived later in this Section. 
Notice that the minimal center of the Clifford group, which is 
required by the group multiplication law, is isomorphic to Zg 
and is generated by e'7l as shown in Ref.^5.. The extension 
of the projective Clifford group by the center 2,% defines the 
minimal unprojected Clifford group which is also finite. 

For U (z^n^^ will denote the corresponding equivalence 
class in by [[/] . A frequently used discrete set of Clifford 
unitaries, whose images [.] generate 1^%^ if we allow them to 
act on arbitrary qubits is 



1 
i 



1 

V2 



1 1 

1 -1 



and (13) 



CNOT= |0)(0|®l2 + |l)(l|®c7i. 

It is obvious that the Pauli group is a subgroup of the Clifford 
group because 



\ c '^„. 



(14) 



The equivalence classes of the elements of which differ 
only by a phase, form a factor-group of ,^„, which we call the 
projective Pauli group, 

and it should play an important role in what follows. Because 
of Eq. (fT4l) the projective Pauli group is also a normal sub- 
group of the projective Clifford group. 



The main results obtained in this paper are based on the 
isomorphism between the projective Clifford group, factor- 
ized by its projective Pauli subgroup [3^,,], and the symplectic 
group Sp2n{'2) over the finite field F2, i.e.. 



^%/[^„]~Sp2„{2). 



(15) 



Although this is a well known result we would like to explain 
it in more detail as it is of central importance for estimating the 
computational power of our topological quantum computer 
The group 5'p2n(2) is defined as the set of invertible 2« x 2« 
matrices with entries from the finite field F2 satisfying the re- 
lation 

SeSp2n{2) ^ S^-M-S^M (16) 
where M is the 2n x 2« block-diagonal matrix 

In order to understand the isomorphism dTSI ) we first note that 
the Pauli operators a, "commute" up to phase factors. This 
means that, if we ignore these phase factors and keep track 
only of the multiplication rules, we obtain an additive group 
with a "modulo-2" arithmetic, which is isomorphic to Zj. In- 
deed, if we identify the Pauli matrices with the vectors 



"or 




'or 


-1 


= I„<8) 


-1 



(0,0) 
(0,1) 



; (73 



(1,0) 
(1,1) 



we could see that the Pauli multiplication rules are trans- 
lated into vector summation over the finite field F2. In other 
words, the one-qubit projective Pauli group is isomorphic to 
the Abelian group 

[^l]~(Z2)2. 

This idea extends naturally to the «-qubit Pauli group 
where we identify the «-qubit Pauli gates with vectors in 
(Fi)^", i.e., the vector space of the 2« -component vectors over 
the field F2, as illustrated by the example 



ail 



)/C73^(l,0,...,l,l). 



Therefore we conclude that the projective n-qubit Pauli group 
is isomorphic to the Abelian group Z?" 



;Z 



2n 



(17) 



Taking into account that the center of the n-qubit Pauli group 
is Z4, generated by the powers of ;, we could compute from 
Eq. dTTl ) the order of the complete Pauli group to be 



■)2n+2 



(18) 



Recall that the Clifford-group elements are structure pre- 
serving, invertible mappings (automorphisms) of the Pauli 
group, i.e. they map ^„ to itself and respect products of Pauli 
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operators. In more detail, if J7 G ^„ and we denote by p the 
multi-index of the Pauli operator cj,, G S^n then 



«-qubit case where p = (pi , . . . ^p^^), q={qi,.. .,q2„) 
and the inner product is given by 



(19) 



where fu{p) is the power exponent of the phase, and we have 
denoted the mapping of the Pauli indices, corresponding to the 
action of the Clifford operator U by 



Su ■■ ¥f 



■F: 



2n 
2 1 



■Su{p)- 



(20) 



While the index mapping ( |20l i essentially expresses the ac- 
tion of the Clifford gate U inside of the projective Pauli group 
[^„], the pair {fij,Su) in Eq. iT% describes unambiguously 
this action on the complete Pauli group The mapping 
( |20l i is furthermore an automorphism of Fj" which could be 
characterized by the following properties: 

1 . Su is linear: 

It follows from the definition, Eq. ( fT9] l. that 



= U*anU , 



where the brackets mean taking the projective equiva- 
lence class. Therefore we can write, using the sign © 
for "modulo-2" addition of vectors. 



U*apec,U 



U*a„a,U 



U*a,,U 



_^Su{p)<BSu{q) 

from which we could derive the relation 

Su{p®q)^Su {p) © Su {qj , V£, £ e ¥l". 
This means that Su £ GL2„{2). 

2. Su is symplectic: 

Working with the projective Pauli group [i3^„] = Fj" 
is convenient, however, we loose the information 
about commutation or anti-commutation of the Pauli- 
operators, e.g. it is impossible to distinguish the two 
cases 



i (T3 • 02 



Ol. 



Nevertheless, it is possible to keep track of the corre- 
sponding commutation relations. If we define the inner 
product of two vectors p = {pi,p2), q ~ {q\,q2) G ¥\ 
to he p * q ~ p2 ■ q\ it is easy to check that 



(-1) 



p*q 



^p®q> 



where C7„ denotes the Pauli-operator associated with a G 
1- 



F?. This generalizes in a straightforward way to the 



i=l 

Whether two operators Op^O^j G ^„ commute or not 
can be inferred from the symplectic form 

« [p^g) ^ p*q-q*p = p^ -M-q, 

which could only be equal to or 1, where M is given 
in Eq. ( fT6l ) and p^ is the transpose of p. It is easy to see 
that ~ ~ 

ap-CT, = (-l)"(M)(j,.(j^, (21) 

from where we can conclude if the two operators Op and 
Oq commute {(o{p,q) ~ 0) or anti-commute {(o{p,q) — 
1). 

Now we can prove that if the operator U in Eq. ( fT9l l be- 
longs to the Clifford group then Su must be symplectic, 
i.e.,yp,qe¥f 



C0{p,q) ^C0{Su{p),Su {q)). 



(22) 



Indeed, let us multiply Eq. ( |2TI ) by U* from the left and 
by U from the right. We have 

U*OpUU*o,,U = i-l)'^^P-3-)u*OgUU*OpU. 

Using Eq. (T% and canceling the phase factor 
lfu(p)+fu{q) on both sides we obtain 

Cysuip)<ysa(q) = i-^r^^'^^<ySaiq)^Suip)- (23) 

On the other side, however, if U G'^n as we assumed, 
then both Osij(p) and Os^[q) are Pauli operators as well, 
so they should also satisfy Eq. ( l2Tl i. i.e.. 



,ip)^s,iq) = (-l)'»(^''(£)'^''(2))cT,,(,)CT,,(rt. (24) 



Comparing Eq. ( [23] ) with Eq. ( l24l i we derive Eq. ( |22] | 
which implies that Su is symplectic, i.e., 

Sl-M-Su=M => SueSp2„{2). 

3. Tlie kernel of U ^ Su is the projective Pauli group 

Consider the mapping between the Clifford and the 
symplectic groups 

S:%^Sp2„i2),U ^Su, 

which is a homomorphism between the two groups %j 
and Sp2n{2) both with matrix multiplication as group 
concatenation. The kernel of this homomorphism, i.e. 
the group of elements in which are mapped to the 
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identity in 5/?2n(2) is generated by those unitaries U 
which satisfy 

t/* (J;i (gi • • • a^„ t/ = ± «) • • • «) C7;„ (25) 

for all choices of indices ji, ■■■,]„ G {0,1,2,3}. Be- 
cause of Eq. (O all Pauli operators satisfy this equa- 
tion, hence 

[^„] cKer5. 

On the other hand, since the Pauli operators form a basis 
for the linear operators on we can express any 

unitary C/ as a linear combination 

Now suppose that there is a [/ satisfying Eq. ( IZSl l with 
more than one terms in the sum above. Then it is easy to 
see that there is a Pauli operator C7y such that the expan- 
sion of U*(JyU consists of more than one Pauli opera- 
tors which contradicts to the assumption that dZSl ) is sat- 
isfied. Therefore the projected Pauli group is not only a 
subgroup of KerS, but also 

Ker5c[^„] ^ Ker5=[^„]. 

Given that the image of the mapping S : U ^ Su is 5'p2«(2) 
and that Ker5 = [^„] we conclude that 

~ ^^„/[^„] ~ Sp2„i2) . 

This important fact allows us to use the results for the order of 
the finite symplectic group 

|5p2„(2)|=2"'n(4^'-l), (26) 

which is proven in Appendix |A] to derive the order of the 
projective Clifford group. It follows from Eq. dTsl l. and the 
fact that [^n] has order 2^", that the order of the projective 
Clifford group is 

l^-^,,! =2"'+2"P](4^'-l). (27) 

This result will be used in Sect. |VT] to estimate the compu- 
tational power of the Ising-anyon topological quantum com- 
puter. 

V. BRAIDING GATES AS CLIFFORD GATES: 
MONODROMY ANALYSIS 

A crucial observation for the results obtained in this paper, 
as well as for topological quantum computation in general, is 
that the Pauli group ^„ for n qubits realized by 2n + 2 Ising 
anyons coincides with the representation of the monodromy 
subgroup of the braid group 

EE Image (.^2n+2) 



Because the monodromy group is in general a normal sub- 
group of the braid group it follows that all braiding gates, 
which by definition leave the monodromy (hence the Pauli 
group) invariant, are actually Clifford gates, i.e.. 

Image (^2n+2) C 

In order to give a more comprehensive explanation of this im- 
portant result we shall recall the standard algebraic definition 
of the monodromy subgroup in terms of the exact sequence of 
group homomorphisms 

where is the monodromy group (also called the pure braid 
group), ^„ is the full braid group and is the symmetric 
group. In more detail, recall that a braid can be described ge- 
ometrically as a bunch of strands connecting the initial config- 
uration of some indistinguishable particles (arranged for sim- 
plicity along a line) with the final configuration. If we forget 
about the strands and look at the final positions of the parti- 
cles, arranged along the same line, then the final configuration 
will be simply a permutation of the original particles. As an 
illustration we show in Fig.|4]the action of the projection map 
on the two-qubit CNOT gate of Refs. 8,9, which is an element 
of the Ising-model representation of Therefore we can 




FIG. 4: Example of the projection mapping K acting on the spe- 
cific braid element CNOT= R:j^RaR3RiR5R4R^^ G ^e, where 
R, = Rij+i in the notation of Ref. |3, producing the permutation 
(216543) e J?6 . 

define the natural "onto"- mapping 

which could be easily seen to be a group-homomorphism. 
Next, as is well known, the kernel of any group homomor- 
phism is a normal subgroup and the factor-group is isomor- 
phic to the symmetric group 

^„/Ker7i:~^„. 

Therefore we may define the monodromy group as the kernel 
of the projection n, i.e., 

= KerTT. 

Geometrically this definition means that a monodromy trans- 
formation is any braid transformation which does not permute 
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the final positions of the particles compared to the initial ones. 
Pictorially this could be described by such braid transforma- 
tions in which one or more particles are transported along 
closed loops around one or more other particles (so that each 
particle arrives at the end at its original position) as shown in 
Fig. |5] It is not difficult to see that these pure braid trans- 



M 




FIG. 5: Decomposition of a monodromy transformation in which the 
particle with coordinate 775 makes a complete loop around particles 
with coordinates 772, r]3 and T74, in terms of the monodromy genera- 
tors Aij. 

formations can be generated by the elementary monodromy 
transformations in which one particle, say with index j, is 
transported along a complete loop around another particle, 
say with index /, where ; < j. Thus the elementary genera- 
tors of the monodromy group could be given the following 
presentation^^ 

Aij = RjliRjlx ■ ■■ ^7+2^7+1 Ri+iRi+2 ■ --Rj-iRj-i 

= U^'R^Uij, where C/;, = J] (28) 

k=i+l 

and 1 < / < y < «. In other words the monodromy group 
could be considered as the smallest invariant subgroup of 3§n 
which is generated by the squares of the elementary braid 
generators^ Rj, /= 1 ,...,«— 1 . As an illustration of Eq. (l28l l 
we display in Fig. |6]the monodromy transformation in which 
the particle with coordinate 775 is transported along a com- 
plete counter-clockwise loop around the particle with coordi- 
nate ri2- 

Alternatively, the monodromy group could be abstractly de- 
fined by the generators Ajj, 1 < / < j <n satisfying certain re- 
lations, however we will skip that as we shall only need here 
the expressions ( |28] ) in terms of the braid generators. 

It can be seen that in the braid-group representation of the 
Ising model the Pauli group is represented in terms of mon- 
odromy operators. Indeed, the (73 gate acting on the ;-th qubit 




TJ, rjz 1, % 76 




'^25 -^45^34 ^23 ^34^^ 



FIG. 6: The monodromy generator A25, representing the elemen- 
tary monodromy transformation in which the particle with coordinate 
r]5 is transported in a counter-clockwise direction along a complete 
loop around the particle with coordinate r]2. expressed in terms of 
the braid generators. 



is simply 

■^') ^ = }ri®_-®h.®a^ (29) 

/— 1 n—i 

while the squares of the non-diagonal braid generators give 
^ = I2 ■ ■ ■ II2 '^g2 I2 , (30) 

i—\ n—i— I 

where 1 < / < « — 1 . The last two squares are more special: 




n 



It is obvious from the above equations that all «-qubit Pauli 
gates can be expressed (up to phases which are powers of /) in 
terms of products of squares of the elementary braid genera- 
tors. For example, the C72 gate acting on the last qubit is 

and the 02 gate acting on qubit (« — j) can be expressed as 
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The above equations also suggest that we might need some ad- 
ditional phase factors in order to construct all elements of ^„ . 
To clarify this point we will construct the element il explic- 
itly in the positive-parity representation (a similar construc- 
tion can also be given for the negative-parity representation 
because both representations are computationally equivalent, 
as shown in Ref. ITll) . More precisely, we will prove below 
that 



R 



/n(n+l.+) 



2n 



^7, 



R 



(n+l,+) 
2n 



{4: 



2«+l 



(33) 



The left hand side corresponds to a composition of two braids, 
where particle 2n moves along a complete loop around par- 
ticles 2n + 1 and 2« + 2 and likewise particle 2n + 1 moves 
around 2« + 2. Therefore Eq. ( l33T l represents an element of 
the monodromy subgroup and we conclude that 



{±I,±/I} C Image 



2n+2) 



To verify ( [33] l we first observe that the last two unprojected 



braid matrices /Jj", and R2I+1 satisfy, according to Eq. (|6]l 
and the anticommutation relations of the 7 matrices, the fol- 
lowing relation 



1) 

2n+l 



(34) 



Indeed, using the properties of the y matrices, it is easy to 
see that the squares of the unprojected braid matrices Q are 
simply 



■ in+l) (n+l) 

-'7) 



1 < y < n - 1- 



(35) 



Applying Eq. ( |35] | for j = 2« + 1, and using the anticommu- 
tation relations (|5]l for the 7 matrices, we have for the LHS of 
Eq. (m 



LHS = 2 ^ 72n+ 1 ) ( - 72»+ 1 72»+2 ) (I - 72n 72«+ 1 ) 
= '(-')72n+l72n+2 = ' (Rln+l)^ ■ 

The fourth power of the unprojected braid matrices is always 
I because applying twice Eq. dSST l we have 



R 



(n+l 



Jn+l) (n+l) (n+l) (n+l) 

1} 0+1 0+1 ' 



which also implies that the fourth power of the projected braid 
generators (fTOl i is I. Finally, the validity of Eq. ( |33] | follows 
from the fact that the projectors (|9]l commute with the group 
multiplication law (see Eq. ( [38] l below). 

Now, because the squares (^^"^^'^^^ of the elementary 

braid generators, as well as the element ;I, belong to the 
monodromy group representation, we conclude that the Pauli 
group is naturally represented in the Ising model TQC as 
a subgroup of the monodromy group for 2« + 2 Ising anyons, 
i.e. , 



On the other hand, as we will show below, the generators ( 
of the monodromy group could be entirely expressed in terms 
of the Pauli operations for n Ising qubits so that 



Image 



2«+2 



) C ^„ ^ S^n = Image 



2«+2 



). (37) 



In order to prove the group isomorphism dJTj l we will show 
that the monodromy generators of the Ising-model repre- 
sentation of ^2n+2 could be expressed (up to powers of /) 
in terms of the squares of the elementary braiding matrices 

, which, according to Eqs. (EB and 

( |32] i. are completely expressible in terms of Pauli operators. 

Before we prove this let us emphasize once again one im- 
portant detail of the Ising representation of the braid group. 
The projectors (|9|l commute with the group multiplication law, 
i.e.. 



^(n+l.±)^("+l.±) ^j^(n+l)j^(n+l)p[n+i:) ^^g^ 



SO that, if we want, we could work all the time with the unpro- 
jected matrices and project them at the end. 

Using Eq. ( l35T l as a base we could prove by induction that 
the unprojected monodromy generators are 

Aw = /(-l)'"*7i"^''7/"^^', 1 < < Z < 2« + 2. (39) 



Indeed, assume that Eq. ( |39] l is fulfilled for given k and I where 
/ > k. Then 

Ajt,/+i = Ri ^AkiRi 

= ^(I + 7/7/+i)('(-1)'"'7aT/) (I-7/7/+1), 

where we have used that the inverse of the unprojected braid 
generator Rj is Ry ' = e""^/"* (I + 7/7/+i ) / V2. Then using the 
anticommutation relations (|5]), to simplify expressions such as 

7/ 7/+ 1 7*: 7 7 7/+ 1 = 77+1 7t 7+1 = -77^, 

we obtain A/^j^i ^ — l)'"'^7t7+i which completes the 
proof of Eq. ( [39] l. Next, plugging between jk and 7/ the 
identities Yk+iYk+i, Yk+iYk+i until 7_i7/_i and re-expressing 
7^7*:+! '{Rk)^ we can write the unprojected monodromy 
generators as 

Aki = -(-/)'-'•+' [Rkf- {Rk+i? ■ ■ ■ {Ri-if [Ri-if . 

Finally, using the fact that projectors commute with the group 
multiplication, as in Eq. ( |38] ), we obtain e.g., the monodromy 
generators in the positive-parity representations (and similar 
for negative parity) as 



i-k+i 



where 1 < k < I < 2n 



R 



k+l) 



W-2)'W-i)'. (40) 



2 and Rt ^ Rl"^'^Pi"+'^ 



Going 

back to Eqs. ( |29] l, ( [30l ), ( [3T] i and ( t32] i we conclude that be- 
cause {RtY G for all k and the element il E S^n as well 
(see Eq. |33]l). the generators of the monodromy representa- 
tion certainly belong to the Pauli group so that 



C Image (^2n+2) • 



(36) 



At, e 



Image (.^2«+2) C 



(41) 
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Combining Eq. ( l4TT i with Eq. (l36T l we finally prove that the 
Ising-model monodromy group representation of ^2n+2 
completely equivalent to the n-qubit Pauli group for the 
Ising-model topological quantum computer. 

Remark 1: It follows from Eq. 071 ) that the center of 
the monodromy group representation coincides with the cen- 
ter of the n-qubit Pauli group which is generated by il, 
i.e.. Center {Image{^2n+2)) = Center{^„) = Z4. Combin- 
ing this with the algebraic fact^ that the permutation group 
■^in+l is centerless for n > I, so that the center of the braid 
group coincides with the center of its monodromy subgroup, 
and with the fact that S^2n+2 is represented faithfully in the 
Ising-model representation of ^2n+2 for n>2, as proven in 
Appendix\C\ we conclude that 

Center (Image (^2«+2)) = Center (Image (.^2n+2)) 

= {±I,±iI}, n>2. (42) 

Remark 2: Because the monodromy group ^2n+2 
Ising-model representation of the braid group S§2n+2 exactly 
coincides with the n-qubit Pauli group this automatically 
means that the image of the braid group is a subgroup of the 
n-qubit Clifford group '^S'n 

Image (^2,1+2) C "^n, 

i.e., all quantum gates that could be implemented by braiding 
in the Ising-model TQC are in fact Clifford gates that stabilize 
the Pauli group. Unfortunately, the converse is not true, i.e., 
it appears that not all Clifford gates could be implemented by 
braiding Ising anyons. 

Furthermore, we could use the group isomorphism dJTl l be- 
tween the monodromy group Image (./#2«+2) and the Pauli 
group to derive the order of the image of the braid group 
■^2n+2- Recall that in general the factor-group .'3§2n+2/ •■^2n+2 
is isomorphic to the permutation group ,S^2n+2- In our case it 
appears that for « > 2 this permutation group is represented 
faithfully (see App. |C] for a proof) so that the order of the 
Ising-model braid-group representation is 

|Image(^2n+2)| = |Image(^2«+2)| |Image(^2n+2)| 

= 22"+2(2«+2)!, «>2, (43) 

where we have used that because of Eq. ( |37] ) the order of the 
monodromy group is equal to the order of the full «-qubit 
Pauli group ,?3^„ (including the center, cf. Eq. (|42]|). which is 
22;i+2 ^jjj j-jj^j. jjjg order of the faithful representation of S^2n+2 
is (2« + 2)!. For n = 1 the order of the image of the mon- 
odromy group is again given by the order of the Pauli group 
22«+2 = 15^ however, the representation of the permutation 
group o5^4 is not faithful. Instead, the third generator is equiv- 
alent to the first one, i.e., 7?^+' = 7?'+' or R^^^ = i (^i"') 
and therefore the factor-group gives a faithful representation 
of the permutation group ^3. Thus the image of the braid 
group for 4 Ising anyons is 

|Image(^4)| = |Image(^4)| |Image(^3)| 
= 16 -3! = 96, 



which exactly coincides with the result of Ref . where the or- 
der has been directly computed using the Dimino's algorithm. 
The result, 2^6! = 46080, for the order of the image of the 
braid group computed from Eq. ( l43T l also exactly matches 
that obtained by the Dimino's algorithm in Ref.[^ 

Finding the order of the image of the braid group ^2n+2 is 
a central result in this paper which is similar to the results of 
Ref. 0, where the order of the image of the braid group de- 
pends on the parity of «, yet it is different. The reason for the 
independence of the order of the representation on the parity 
of n in our case is the presence of the additional phase fac- 
tor e'7 in Eq. which certainly changes the center of the 
representation of the braid group. 

VI. THE PROJECTIVE CLIFFORD GROUP AND 
THE ISING REPRESENTATION OF THE BRAID GROUP 

In this Section we shall prove that not all Clifford-group 
gates could actually be implemented by braiding of Ising 
anyons. To this end we shall demonstrate that the order of the 
projective Clifford group, i.e., the Clifford group factorized 
by its center, is much bigger and grows much faster with n 
than the order of the image of the braid group. Because global 
phases are irrelevant for quantum computation it makes sense 
to compare the orders of the two finite groups factorized by 
their centers. 

Let us first compare the single-qubit braid group S§4 and 
the single-qubit Clifford group The order of the Ising- 

model representation of ^4 is 96, however if we factorize it by 
its Z4 center we have |Image(,^4) /Z4I = 24. This precisely 
coincides with the order of the projective Clifford group for 
one qubit obtained from Eq. ( |27] |. 

For two qubits the orders of the representation of the braid 
group and its projective version are 

|Image(i^6)| = 46080 ^ \lmiig&{^e)IM = 11520 

and again the second number precisely coincides with the or- 
der of the two-qubit projective Clifford group obtained from 
Eq. ^ 

l^-Tzl =22'+4(4- 1)(42- 1) = 11520. 

While the coincidence of the orders of the one-qubit braid 
group representation and the one-qubit Clifford group seems 
natural, the one for the two-qubit groups is more sub- 
tle. It could be understood by the occasional group 
isomorphis m^*^'^^ 

and this remarkable mathematical fact could eventually ex- 
plain why the entire Clifford group for two qubits can indeed 
be realized by braiding, while for any n > 3 the Clifford group 
is much bigger, and grows much faster with «, than the image 
of the braid group. Indeed, the orders of the image of the 
braid group (factored by its center) and of the projective Clif- 
ford group for n = 1 , ... 5 number of qubits are compared in 
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Table J] Already for 3 qubits the order of the projective Clif- 
ford group is about 40 times bigger than the order of the image 
of the braid group factored by its center. For n —4 and 5 the 

TABLE I: Comparing the orders of the image of the braid group 
^2n+2' factorized by its center Z4, and of the projective Clifford 
group .^"rfn. 



# of qubits n 


1 


2 


3 


4 


5 


|Im(^2«+2)/Z4| 


24 


11520 


2580480 


oc 0.9 X 10'' 


o<0.5x 10'^ 




24 


11520 


92897280 


oc 1.2 X 10'^ 


CK2.5X 10''' 



discrepancy is even more dramatic: the order of the projective 
Clifford group is increasing much faster with n than the order 
of the image of the braid group. In Fig. [7] we have plotted 
the logarithm of the ratio between the order of the projective 
Clifford group and the order of the braid group for number of 
qubits up to 100. To understand more intuitively Fig.|2]notice 




10 20 30 40 50 60 70 80 90 100 
number of qubits n 

FIG. 7: Comparing the order of the «-qubit Clifford group %, with 
that of the image of the n-qubit braid group .^2«+2- The logarithm 
of the ratio of the two orders is plotted on the vertical axis versus the 
number of qubits n. 



of Ising anyons. These are the embeddings of the two-qubit 
SWAP gate into the n-qubit system. Indeed, if we have them 
then we could construct by braiding all embeddings of the CZ 
gates because by physically exchanging pairs of Ising anyons, 
representing the qubits, we readily obtain CZ.SWAP as illus- 
trated in Eq. ( |45T l below for the two-qubit case. In general, by 
exchanging the j-th and {j + 1 )-th pair of anyons, with coordi- 
nates (J72j- 1,772;) and (J72j+i, J72;+2) respectively, we obtain 
(cf. g^) 

where CZ'^^ -'^'^ is the Control -Z gate acting on the j-th qubit 
(the control qubit) and {j + l)-th qubit (the target qubit) and 
SWAP(^'^+') is the two-qubit SWAP gate exchanging the ;-th 
and {j + 1 )-th qubits. The Phase gate P in ( fT3] l acting on qubit 

j is given by /?2/-i^'- Moreover, if we have all the SWAPs 
then we can construct Hadamard gate acting on arbitrary qubit 
by swapping the Hadamard gate 



(«) - 



W"i = R 



■p(«+l. + )n(n+l. + ) 
^2«+l ^2n 



(4: 



+1,+) 

2n+l 



acting on the last qubit to the qubit where it is needed. 

Thus we see that if we could construct all embeddings of 
the two-qubit SWAP gate then we could construct the entire 
Clifford group for n qubits. 



VII. SOLUTION FOR TWO QUBITS: BRAIDING 
IMPLEMENTATION OF THE SWAP GATE 

The problem with the implementation of all Clifford gates 
could be most easily illustrated on the example of the two- 
qubit SWAP gate^'\ which is a Clifford gate that can be ex- 
pressed in terms of CNOT and Hadamard gates ^ as shown in 
Fig. [H Imagine a two-qubit system constructed from 6 Ising 



that the leading term in the order of the projective Clifford 
group is 

I ^"^^ I ~ 2"'+2" fl oc 22"'+3«, 

j=i 

while (using the Stirling formula) the order of the image of 
the (projective) braid group is 

|Image(<^2«+2)| =22«(2« + 2)!oc 22"+(2«+2)iog2(2«+2) 
and therefore the leading term in the ratio is 

|^5^"^„|/|lmagei^2«+2| - 22"'+«-(2«+2)iog2(2«+2) ^ 22"', 

so that the logarithm of this ratio still grows quadratically with 
«. 

It is not difficult to see which gates are missing, i.e., which 
are the Clifford gates that cannot be implemented by braiding 





H 












H 









H 



FIG. 8: SWAP gate in terms of 3 CNOT and 4 Hadamard gates 

anyons, with coordinates 771 , . . . , 775 in the plane, in such a way 
that the first pair (with coordinates 771 and 772) forms the first 
qubit, the second pair (with coordinates 773 and 774) - the sec- 
ond qubit and the last two qubits with coordinates 775 and 776 
form the inert pair as shown in Fig.|9] Without restriction of 
generality we could choose to work in the spinor representa- 
tion S+ with positive total parity because, as shown in Ref.— , 
both inequivalent spinor representations S± of the braid group 
■^2n+2 are computationally equivalent. The two-qubit compu- 
tational basis in the spinor representation S+ can be written in 
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FIG. 9: Two-qubits configurations in terms of 6 Ising quasiholes cor- 
responding to the positive-parity representation of the braid group 



terms of Ising-model correlation functions as 

|00) = {a+a+ a+a+ a+a+) 
|01) = (cT+(T+ a+a (T+cT-) 

1 10) = {a+G- a+a+ a+G-) 

111) = (c7+(7_ a+a (7+C7+) 



(44) 



Now, let us try to swap the two qubits by simply exchanging 
the two pairs comprising the qubits. This could be done by 
the following braid transformation (see Ref.— for the explicit 
form of the braid matrices in this basis) 



.,(2) _ p(3,+)p(3,+) (3,+) (3,+) 
^#l«#2-^2 ^3 ^1 ^2 



10 
10 

10 

-1 



(45) 



The resulting unitary transformation has one extra minus on 
the last row as compared to the standard SWAP gate' (we skip 
the global phase /). This is, however, natural because this mi- 
nus sign appears in front of the state 1 1 1 ) that is the state in 
which both Ising anyon pairs at (771 , 772) and at (773, T74) share 
Majorana fermions and therefore exchanging the two qubits 
is equivalent (especially after fusing 771 ^ 772 and 773 774) 
to the exchange of two fermions which certainly produces a 
minus sign. This raises the reasonable question whether the 
SWAP gates is at all realizable by braiding. 

On the other hand, we could have chosen a different basis 
of Ising (7 fields to represent two qubits, e.g., the first pair 
(with coordinates 771 and 772) forms the second qubit, while 
the second pair (with coordinates 773 and 774) forms the first 
qubit and the last two qubits with coordinates 775 and 775 form 
the inert pair again as shown in Fig. \TU\ Then the two-qubit 




correlation functions as 

100) ' ^ 

101) ' ^ 

1 10) ' ^ 

111) ' ^ 



,C7+C7+ (T+C7+ C7+CJ+) 
'C7+C7_ <y+<y+ C+CT-) 
'C7+C7+ <J+<y C7+(7_) 



(46) 



Obviously the basis ( |46] | could be obtained from the basis (l44l i 
by the action of the transformation matrix 



1(2) 



10 
10 
10 
1 



; SWAP. 



(47) 



^3 ^4 



which simply means that the bases (l44l) and ( l46b describe 
equivalent representations of the braid group ^g. 

The question whether the SWAP gate is implementable by 
braiding is equivalent to that if the product 

f/#i«#2A(2)=/diag(l, 1,1,-1), 

which is (up to phase) equal to the Controlled-Z gate, is ex- 
pressible in terms of the braid-group generators. Indeed, using 
the explicit form of the generators from Ref.[Tll we could eas- 
ily check that 

/^f^+)7^f^+)(/^f+')"^=diag(l,l,l,-l)=CZ, 

which constructively proves that the two-qubit CZ and SWAP 
gates (hence, all Clifford gates for two qubits) are efficiently 
realizable by braiding of Ising anyons. 



VIII. CONCLUSION 

In this paper we have demonstrated that the Pauli group in 
the Ising-anyon topological quantum computer exactly coin- 
cides with the image of the monodromy group. Therefore, all 
quantum gates that can be implemented by braiding of Ising 
anyons are Clifford gates that stabilize the Pauli group. An- 
alyzing the structure of the Clifford group for « -qubits and 
comparing to the structure of the braid-group representation of 
the Ising model, we conclude that not all Clifford gates could 
in general be realized by braiding in this model. The gates that 
are missing are the embeddings of the two-qubit SWAP gate 
into a n-qubit Ising systems. Thus, not only the Ising-anyon 
topological quantum computer is not universal for topologi- 
cal quantum computation but even the entire Clifford group is 
not implementable by braiding. Nevertheless, the braid-group 
representation of the Ising TQC contains all of the entangling 
gates that could be used to analyze topological entanglement 
and topological protection of qubits and quantum operations. 



FIG. 10: Two-qubits configurations conesponding to the positive Acknowledgments 
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APPENDIX A: ORDER OF THE SYMPLECTIC GROUP 

Following Refi^^ we could find the order of the symplec- 
tic group Sp2n{q) over the field F^, by counting the number 
of different ways to chose a symplectic basis, that we shall 
order as {e\ ,.. . ,e„;/i , . . . ,/„}, of the vector space F^", i.e., 
the bilinear form in this basis is (e,-,/i) = 1 = — (/,-,e,) and 
{e.Jj) =0= {fj,ei) for i^j. 

1 . Choose ei : Any non-zero vector in F^" could be chosen 
as ei. The number of vectors in F^" is q^" so that the 



number of non-zero vectors is #e 



1. 



2. Choose /i: The vector /i might be chosen among the 
vectors which are not orthogonal to ei. The dimension 
of the subspace {v e F^"|(v,ei) = 0} is 2n - 1 (at this 
point we could think of ei as being an element of an 
orthonormal basis of F^") so that the number of the or- 
thogonal vectors is q^"^ ' . 

3. The number of vectors which are not orthogonal to ei 
is thus q^" — q^"^^ = (<? ^ 1)<7^"^' and all they come in 
scalar multiples such as, fi , 2/i , . . . , — so that 
the number of different choices for fi is #/i = q^"^^. 

4. Therefore the number of different choices of ei ,/i is 

#{eufi}^iq'"-l)q"'-'- 

5. The rest of the space, which is spanned by the basis 
vectors {e2, ■ ■ ■ ,e„;/2, . . . ,/„} is the vector space F^"^^ 
of dimension 2 (n — 1) so we can apply again the above 
result to obtain 

(after {ei,fi} have been chosen). 

6. Induction: continuing in this way we finally obtain 

\Sp2.{q)\ - n {q'-'-l)q^'-'=q"'U{q^'-l)- (Al) 
;="i i=i 

Applying Eq. (lAll i for q = 2we obtain Eq. 



APPENDIX B: ISING BRAIDING GATES AS SYMPLECTIC 
TRANSFORMATIONS 

We will determine the symplectic matrices corresponding 
to the Ising model exchange matrices via relation (flSl l. To 
this end we first phrase the Ising model quantum computer 
abstractly in three steps: 



• Starting point is the representation of the Clifford alge- 
bra in terms of the Ty"^''- According to ^ we define 
operators Rj which act on where = C^. 

• The subspace of jf^J"+^ of positive parity, i.e. the lin- 
ear span of the set 



{\x-i...x„+i) : xit 



)a:„+i =0}, 



will be denoted by J^f "+'. We will identify 
with Jf^'" via the correspondence 

\xi...x„z) \xi...x,j) 

where z is chosen such that the left hand side is an ele- 
ment of i.e. z=xi®---S)x„. 

• Finally, the action of the Rf on \xi . . .x„) G is ob- 
tained by applying Rj to \x\ . . .Xnz) and calculating the 
corresponding vector in Note that this step is 

only well-defined, because Rj is reducible with respect 
to the decomposition = © 

The last two steps describe the action of the projection P^'^ in 
section Hin 

Now, instead of focusing on the correspondence between 
basis vectors, we will concentrate on the correspondence be- 
tween the Pauli groups associated with and 
The operators 



(Bl) 



if J 7^ 3 and 



(B2) 



where we denote by |"+ the restriction to Jif®"^^, are essen- 
tially the Pauli operators on This follows from the fact 
that these operators are also reducible with respect to the de- 
composition = ® ^®"+' which allows us to 
calculate their action on ,J^®" just like we do for the Rj . 

Since both, Rj and R^ are Clifford unitaries in their respec- 

tive domains, it is easy to determine the action of RJ on Tj 
according to the following commutative diagram 



{R+)*zfR+ 



R/Gj Gg Ri 



T 

n 



=1 



where we have used the short hand notation 



Note that the C7g"^'' in the lower left corner depends on the 
af, we have 5 = 1 if ; = 1,2 and 5 = if ; = 0,3. 
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The action of the operators Rj on Pauli operators is 
equivalent^ to the following Clifford unitaries 



1 
i 



and similarly for T, . The relations for can be calculated 
in the same way, but we will skip the proof and give only the 
result. If A: < n — 1 we have 



(k) Rtn^ 



n 



1 



10 0/ 

1 -/ 
0-/10 
/ 1 



These Clifford gates map Pauli operators according to 



(72 



(72 



(k) «2„+l 



The case k = n gives 



n4" 

\r^k,n 



(n) *2n+I (n) 



"1 



(7i (g)ll 
1(E) Gi 



(73 ® (72 
(72® (73 



)(72 ■ 



) (72 



It is straightforward to determine the action of the /?t for all 
j = l,...,2n— 1. Since the corresponding Rj do not effect 
qubit n + 1 we get non-trivial relations only for qubits j and 



R 



2j-l ■ 

Rtr- 



,(;■) 



Jj) 



''3 '2 



JJ) 



rU) 



JJ) 



Jj) 



Jj) 

2 

Jj) 



The action of the last two operators can be calculated by uti- 
lizing 



It follows that 



r^k 



;I 



W Jk) 



n 

r^k 



3 ' 



which can be verified as follows 



(R+ Yt^''^R+ - (R+ Yn^^'n 



- (}?+ 



2)1+1 1 + 

- (^2h+i) ^2/7+1 r+ 
= C7«c7i-^' 



)C72( 



)(7i 



From these equations we can read off the corresponding sym- 
plectic matrices 



52, 



i-1 





Mi-y) 








1 — 





1 







1 












\(n-i) _ 



for 1 < 2/ - 1< 2« - 1 and 

















1 










52,= 





1 


1 1 













1 











1 


1 


1 












l2(n-i-l) 



for 1 < 2/ < 2n — 2. The last two matrices read 



Sin 



Sin 



+ 1 



1 

1 


1 1 
1 1 




1 
1 


1 1 
1 1 


1 

1 




1 
1 










1 1 



1 

1 


1 1 



1 1 
1 1 


1 
1 


1 




1 
1 


1 1 

1 1 


1 

1 


1 1 
1 1 




1 1 
1 1 


1 1 
1 1 


1 

1 
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APPENDIX C: FAITHFULNESS OF THE 
REPRESENTATION OF THE PERMUTATION GROUP 

In this appendix we will use the symplectic description of 
the braid generators given in Appendix |B] to prove that the 
equivalence classes (i.e., the cosets of the monodromy sub- 
group) of the elementary braid generators are all different and 
hence they form a faithful representation of the symmetric 
group. According to section|V]we have the relation 



which changes between both bases. The new symplectic ma- 



Image(^ 



nn+2 



) /Image (.^2n+2) C Sp2n{2) 



and therefore the advantage of the symplectic description is 
that symplectic matrices already describe classes of braiding 
operators with respect to the equivalence relation in which two 
elements are equivalent if they differ by a monodromy trans- 
formation only and therefore it is natural to compare these 
matrices directly. 

Recall that a representation of the permutation group ^2/7+2 
is by definition generated by 2n + 1 matrices 



satisfying the relations 

SjSj+iSj 
SjSk 



l<j <2n + \ 



Sj+iSjSj+u 
SkSj, \k-j\>l, 



and this representation is faithful if and only if each generator 
Sj is independent from the others with respect to the group 
multiplication. We denote the elementary generators of the 
abstract group ^2h+2 by TZj. It is understood that each Sj rep- 
resents Tlj. 

It is easy to see that the matrices Sj, that we derived in 
Appendix |B] indeed satisfy the above relations. To further 
analyze the representation given by the Sj it is convenient to 
change the basis in the symplectic space according to 

(0_^, 1 , 0, 0^_^) ^ (0^__0, 1 , 

2/ 2(H-i-l) 2i 2(H-i-l) 



(0^__0, 0, 1 , 0^_^) ^ (0^__0, 0, 1 ^ 

2i 2(n-i-l) 2i 2(n-i-l) 

To obtain the matrices Sj in this new basis we have to multiply 
them from both sides with the (self-inverse over F2) matrix 



1 
1 












1 1 
1 1 


1 
1 
















1 1 
1 1 


1 1 
1 1 




1 
1 



trices Sj = T ■ Sj ■ T are 



Si 



1-1 



fori <2i-l < 2/7-1, 


















1 









1 












h{n-i) _ 





"l2,-l 








Sii^ 





1 

1 
















for 1 < 2; < 2« — 2 and the last two are 





" 1 








. 


. 


1 ' 







1 











1 










1 








1 


Sin = 











1 





1 













. 


. 1 


1 




. 














1 . 



Sin 



+ 1 



1 

1 


1 1 
1 1 


1 1 
1 1 




1 1 
1 1 


1 

1 


1 1 
1 1 




1 1 
1 1 


1 1 
1 1 


1 

1 













It is now obvious that the Sj for 7 = 1 , . . . , 2« — 1 represent 
elementary transpositions. Moreover, it is clear that Si„ can 
never be constructed solely from the Sj for j — 1 , . . . , 2n — 1 
because of its non-trivial last column. Therefore the operators 
§1,..., Sin generate a faithful representation of ^2;!+i ■ 

If n = 1 the two 2x2 matrices 53 and §1 coincide, hence 
we get a faithful representation of S^j. However, for n > 2 
a simple argument shows that Sin+i is independent from all 
other Sj. Indeed, suppose there is a way of writing Sin+i as a 
product of the Sj with j < 2n 



Sin 



+ 1 



■Sa,., l<a;<2«. 



The right hand side corresponds to a permutation n E ^in+i 
whose expression in elementary group generators Tij reads 



71 — 7t. 



«1 



'■a, 



Alternatively, the element n can always be written as 



n^«- ) ( ) , j3, < 2« - 1 
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Therefore we also have the relation 

Now observe that multiplying S2n with 5/s with j < 2« — 1 
simply exchanges rows and columns in 52„. The number of 



I's in 52,7 is left unchanged by exchanging rows or columns. 
Therefore, it is obvious that S2n+\ cannot be constructed like 
that from Sin, because the number of 1 's of both matrices dif- 
fers. Hence, we have shown that the representation of the per- 
mutation group =5^2,1+2 is faithful for all n>2. 
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The operators Rj and P, R differ only by a global phase factor, 
which means that the corresponding elements in the projective 



Clifford group ,0^%, coincide. 



